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RESUMEN 

Se presentan relaciones de escala entre la velocidad virial (V) y la dispersion 
, de velocidades uni-dimensional central (<x ); el radio gravitacional (R v ) y el radio 

CsJ ■ efectivo (R e ); y la masa total (M) y la masa luminosa (M L ) encontradas en simula- 

q^, ciones de /V-cuerpos de fusiones binarias entre galaxias espirales. Estas relaciones 

1} ' son de la forma V 2 oc ag , i? v cx i?f y M oc M^. Los valores particulares obtenidos 

m : para {a,(3,j} dependen del metodo de ajuste utilizado [mmimos cuadrados ordi- 

' narios (ols) o regresion ortogonal de distancia (odr)], el tipo de perfil supuesto [de 

Vaucouleurs (deV) or Sersic (S)], y el tamano del intervalo radial donde se hace el 
ajuste. Los indices a y 7 resultan ser los mas sensibles al procedimiento de ajuste, 
^ , obteniendose para el OLS un promedio (a) \ s = 1.51 y (7)ois = 0.69, mientras que 

para el ODR un (a) dr = 2.35 y (7)odr = 0.76. El indice (3 depende mas de el perfil 
adoptado, con (f3)dcV = 0.13 y (/3)s = 0.27. Concluimos que los remanentes de fu- 
^\ ' siones formados de manera no disipativa tienen un rompimicnto fuerte de homologia 

, cstructural y cinematica. 
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We present scaling relations between the virial velocity (V) and the one- 
dimensional central velocity dispersion (cr ); the gravitational radius (i? v ) and the 



effective radius (R c ); and the total mass (M) and the luminous mass (Ml) found in 
iV-body simulations of binary mergers of spiral galaxies. These scalings are of the 
form V 2 oc <7q, R v oc R@ and M oc M^. The particlar values obtained for {a,/3,7} 
depend on the method of fitting used [ordinary least-squares (ols) or orthogonal 
distance regression (odr)], the assumed profile [de Vaucouleurs (deV) or Sersic 
, (S)], and the size of the radial interval where the fit is done. The a and 7 indexes 

' turn out more sensitive to the fitting procedure, obtaining for the OLS a mean 

(a)ois = 1.51 and (7) ols = 0.69, while for the odr (a) dr = 2.35 and (7) dr = 0.76. 
The (3 index depends more on the adopted type of profile, with (/3)d c v = 0.13 and 
(P)s = 0.27. We conclude that dissipationless formed remnants of mergers have a 
strong breaking of structural and kinematical homology. 

Key Words: GALAXIES: KINEMATICS AND DYNAMICS, METH- 
ODS: NUMERICAL, TV-BODY SIMULATIONS 



1. INTRODUCTION Ellipticals show a number of regularities among 

their kinematical and structural properties that have 
been recognized in the past, such as the Kormendy 
and Faber-Jackson relation, and the Fundamental 



Toomre's (1977) idea that the merging of spirals 
could lead to an elliptical galaxy has found ground 

evidence, both theoretical and observational (e.g., ,. . , 

inr>o o 1 inno\ l^i 1 Plane of Ellipticals (e.g., Kormendy 1977, laber & 

Barnes 1998, Schweizer 1998), although some ques- T , „„„„ ^ , , \ 

/ n 11 r,r>r>r, • o Jackson 1976, Bernardi et al. 2003a, b). 

tions remain open (e.g., Peebles 2002, Chiosi & Car- ' 

raro 2002). Understanding the physical origin of these rela- 

tions is important since they are intimately related to 
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theoretical scalings relations among physical and ob- 
servational quantities depends on our knowledge, for 
example, of the star formation process, the distribu- 
tion of dark matter in ellipticals, and the kinematics 
of remnants of spiral mergers. Given the complex- 
ities of the problem, we restrict ourselves here to 
obtain scaling relations resulting solely from dissipa- 
tionlcss simulations of mergers of spirals. 

In a related paper (Aceves & Velazquez 2005) 
the accumulated effects of spatial and kinematical 
homology breaking on the determination of a Fun- 
damental Plane (FP) like relation for remnants was 
studied, but no detailed examination of how the dif- 
ferent physical quantities involved scaled with their 
observational counterparts. In this work we ad- 
dress this matter in more detail, and include two 
more merger simulations than in Aceves & Velazquez 
(2005). 

In particular, we consider here only quantities 
involved in the virial theorem and we determine their 
dependences with their observational counterparts: 



V 2 oc <t£, R e cx RP 



M<xM£', 



(1) 



where V is virial velocity, a® the one-dimensional 
central velocity dispersion, R g the gravitational ra- 
dius, i? the effective radius (i.e., that enclosing half 
of the luminous matter), M the total mass of the 
system, and Ml the luminous mass. Homology be- 
tween the physical and observational quantities re- 
quires that a = 2, and (3 = 7 = 1 in expression 

A further motivation for this study stems from 
the fact that behavior of the previous relations bears 
direct impact on the estimate of a dynamical mass 
Mdyn in ellipticals (e. g., Padmanabhan et al. 2004). 
It is common to assume in such studies an homolo- 
gous relation of the form M dyn oc <r 2 R c . However, 
as we will see, merger remnants show an important 
deviation from homology that most probably reflect 
the actual situation in ellipticals. Thus this effect 
of non-homology would need to be taken into con- 
sideration when estimating the dynamical mass of 
ellipticals. However, the study of such a a problem 
is out of the scope of the present work. 

This paper has been organized as follows. In Sec- 
tion 2, we summarize the numerical models of spi- 
ral galaxies used. We describe the initial conditions 
for the encounters as well as we provide some infor- 
mation on the computational aspects on the simula- 
tions carried out. In Section 3 we present the method 
used to determine the physical quantities {V, R B , M} 
and their observationally related ones {<jq, R e , Ml,}. 
Also, the fitting procedures used to obtain the scal- 



ing indexes in are indicated. In Section 4, the 
results obtained for {a,(3,j} in JJJ are presented. 
Finally, in Section 5, we summarize our main con- 
clusions. 

2. GALAXY MODELS AND INITIAL 
CONDITIONS FOR ENCOUNTERS 

2.1. Galaxy Models 

Our spiral galaxy models consist of a spherical 
dark halo and a stellar disk component. The contri- 
bution of a central bulge is not considered here. The 
disk profile has the functional form 



Pd(R, z) 



M d 



47Ti?2 



d z d 



■exp(-R/R d )sech 2 (z/z d ) , (2) 



where R d and z d are the radial and vertical scale- 
lengths of the disk, respectively. The vertical length, 
Zd, is randomly taken from the interval (0.1 — 0.2)i?d; 
where R d is obtained as indicated below. 

The dark halo follows a Navarro, Frenk & White 
(NFW, 1997) profile, modified with an exponential 
cutoff: 



Ph(r) 



M h a h 



47rr(r + r s ) 2 



exp 



^200 



. (3) 



with 



exp(g 2 ) 



^qexp(q 2 )Eifc(q) + \ exp(g 2 )Ei(g 2 ) - 1 



where Erfc(x) is the complimentary error function 
and Ei (a;) the exponential integral. The scale radius 
of the dark matter profile is r s , c = 1/q = t^ooAs is 
the concentration, and Mh is the halo mass; r2oo is 
defined as the radius where the mean interior density 
is 200 times the critical density. 

The properties of the disk are set up satisfying 
the Tully-Fisher relation (Tully & Fisher 1977, Gio- 
vanelly et al. 1997). This is carried out by following 
the study of Shen, Mo & Shu (2002, hereafter SMS) 
and using the disk galaxy formation model of Mo, 
Mao & White (1998, hereafter MMW); from which 
we can obtain R d . 

In the MMW framework five parameters are re- 
quired to obtain the radial scale-length of the disk. 
These are the circular velocity V c at r2oo, the dimen- 
sionless spin parameter A, the concentration c of the 
dark halo, the fraction of disk to halo mass m d , the 
fraction of angular momentum in the disk to that in 
the halo jd- We have followed the procedure outlined 
by SMS to construct our galaxy models, and chosen 
an epoch for the formation of disks at a redshift of 
0=1 (Peebles 1993). 
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We have selected only spirals with circular ve- 
locities in the range from 50 to 300 km s -1 , and 
with a disk stability parameter e m > 0.9; where 
£m=V m (GMd/i?d) -1 ^ 2 an d V m is the maximum rota- 
tion velocity (Efsthatiou, Lake & Negroponte 1982, 
Syer, Mao & Mo 1997). From an ensemble of random 
points constructed according to the scheme of SMS, 
and satisfying the previous conditions, we obtained 
the final properties of the 24 galaxies that take part 
in our 12 binary merger simulations. 

In Tabled the particular values for each galaxy 
model constructed are listed; where Ah and are 
the number of particles used in the halo and disk, 
respectively. The last column lists the pericenter 
radius R p for the encounters, assuming that galax- 
ies are point particles. Finally, Hernquist's method 
(1993) was used to set up the particle velocities in 
our self-consistent models. 

2.2. Encounter Parameters 

A large number of simulations, and computa- 
tional resources, would be required to sample the 
parameter space of binary encounters of disk galax- 
ies in order to address the dynamical effects of, for 
example, the pericenter distance and disk orienta- 
tions on the scaling indexes {a,/3, 7} in (T|. 

We decided instead to sample randomly the 
encounters initial conditions, but considering only 
parabolic encounters. The pericenters R p were cho- 
sen randomly in the range of {5 — 20} kpc; values 
that are typically found in cosmological simulations 
and that tend to favor mergers (e.g., Navarro, Frenk 
& White 1995). The particular values of R p are in- 
dicated in Tabled 

The initial separation between two galaxies is 
25% larger than the sum of their corresponding r2oo 
radii. The spin orientation of each galaxy, relative 
to the orbital plane, is taken also randomly. 

2.3. Computational Issues 

The simulations were done using GADGET, a 
tree-based code (Springel, Yoshida & White 2001), 
and run on a Pentium cluster of 32 processors 
(Velazquez & Aguilar 2003). We chose the soft- 
ening parameter for disk particles £d = 35 pc and 
£h = 350pc for dark particles. 

GADGET uses a spline kernel for the softening, 
so the gravitational interaction between two particles 
is fully Newtonian for separations larger than twice 
the softening parameter (Power et al. 2003). This 
corresponds in practice to the numerical resolution 
of our simulations. 

We evolved in isolation the numerical realizations 
of each galaxy for about 2 Gyr, and no significant 



change was appreciated in their density profiles or 
virial ratio. Each binary merger was followed for a 
total time of about 8 Gyr. At this time the remnants 
had reached a stable virial ratio. The typical time 
of arrival to pericenter is about 1 Gyr. Energy con- 
servation was better than 0.25% in all simulations. 
Each simulation took « 2 weeks of wall clock time 
in our PC cluster. 

The center of each remnant was determined by 
the center-of-mass of the 1% most bounded particles. 
We eliminated any residual bulk motion from the 
remnant before computing their properties. 

3. METHOD 

In this section we describe how the different 
physical and observational quantities were obtained. 
Also, the different fitting procedures used to obtain 
{a, /?, 7} in Q are indicated. 

3.1. Physical Quantities 

The virial relation may be written as V 2 = 
GM/Rg] where V 2 is the 3-dimensional velocity dis- 
persion, M the total mass of the system and R g the 
gravitational radius. We estimate these quantities 
as: V 2 = 2T/M and R g = GM 2 /\W\, where T and 
W are the kinetic and gravitational energy of the 
remnants, respectively. 

The total kinetic energy T and gravitational en- 
ergy W were computed from the usual formulae: 

T=i5>i*, ^ = -G££^, (4) 

Z— 1 % 1<J J 

where ry the separation between particles i-th and 
j-th. The summation is taken only over the bound 
particles of the resulting remnant. 

3.2. Observational Quantities 

In order to obtain R c an observational "proce- 
dure" was followed. We fitted an assumed pro- 
file to the surface density profile of the luminous 
mass, E(i?). We choose for this the i? 1 / 4 -profile 
(de Vaucouleurs 1948) and Sersic i? 1 / ra -profile (Sersic 
1968). These are analytical formulae that are com- 
monly used in observational studies of ellipticals 
(e.g., Caon, Capaccioli & D'Onofrio 1993). 

Sersic profile has the form 

S(i?) =S exp[-6(i?/ J R e ) 1 /"] ; (5) 

where b = b(n). This profile reduces to de Vau- 
couleurs one when the index n = 4; for n = 1 an 
exponential profile is obtained. 



SCALINGS IN MERGERS 



3 



TABLE 1 

PROPERTIES OF INITIAL GALAXIES 









Halo 










Disk 










Mh 


' zUU 


A 


j 


c 


Nh 

n 


Ma 


Rd 




Na 






L 1VJ -0J 


[kncl 










L 1VJ -0J 










M01 


4.02 x 10 11 


104.3 


0.063 


0.041 


7.63 


57126 


2.11 x 10 10 


2.4 


0.39 


12000 


13.2 




1 25 x 10 12 


152.3 


0.056 


0.068 


3.84 


177571 


7.78 x 10 10 


5.7 


0.47 


44267 




M02 


1.39 x 10 11 


73.3 


0.021 


0.052 


8.94 


172403 


4.48 x 10 9 


1.2 


0.12 


57436 


13.3 




6 45 x 10 10 


56.7 


0.028 


0.033 


6.61 


80000 


1 56 x 1 9 


1.3 


0.13 


20000 




M03 


8.65 x 10 10 


62.6 


0.078 


0.033 


11.12 


149138 


3.62 x 10 9 


1.7 


0.20 


42093 


6.9 




4.64 x 10 10 


50.8 


0.031 


0.043 


12.38 


80000 


1.72 x 10 9 


0.7 


0.12 


20000 

-^V/V/V/V/ 




M04 


2.36 x 10 11 


87.4 


0.036 


0.029 


9.17 


80000 


8.64 x 10 9 


1.1 


0.19 


20000 


10.4 




1.75 x 10 11 


79.2 


0.048 


0.054 


5.29 


59428 


6.77 x 10 9 


3.3 


0.58 


15671 




M05 


5.42 x 10 10 


53.5 


0.053 


0.114 


15.95 


74829 


4.82 x 10 9 


1.2 


0.19 


20000 


18.6 




6.02 x 10 10 


55.4 


0.034 


0.077 


12.05 


83124 


3.44 x 10 9 


1.0 


0.14 


14245 




M06 


7.49 x 10 10 


59.6 


0.098 


0.059 


7.91 


91546 


5.11 x 10 9 


2.5 


0.35 


25000 


15.5 




6 68 x 10 10 


57.4 


0.078 


0.043 


13.32 


81742 


2.93 x 10 9 


1.9 


0.23 


14327 




M07 


4.88 x 10 10 


51.7 


0.074 


0.064 


11.02 


73643 


3.31 x 10 9 


1.5 


0.25 


20000 


14.1 




9 48 x 10 10 

-/.IV" ' \ IV/ 


64.5 


0.047 


0.103 


11.79 


143062 


6 98 x 1 9 


1.6 


0.31 


42119 




M08 


9.77 x 10 10 


65.2 


0.032 


0.034 


6.56 


245462 


2.39 x 10 9 


1.7 


0.20 


36000 


12.5 




1.02 x 10 11 


66.1 


0.023 


0.012 


7.80 


188124 


1.31 x 10 9 


0.9 


0.12 


19725 




M09 


8.11 x 10 10 


61.2 


0.099 


0.142 


10.87 


150000 


7.22 x 10 9 


4.6 


0.73 


30000 


8.4 




8.33 x 10 10 


61.8 


0.122 


0.095 


10.01 


154050 


7.15 x 10 9 


3.9 


0.56 


29686 




M10 


1.27 x 10 11 


71.0 


0.109 


0.078 


9.42 


240000 


9.14 x 10 9 


4.0 


0.42 


60000 


7.9 




8.76 x 10 10 


62.8 


0.088 


0.089 


9.23 


165992 


7.42 x 10 9 


2.9 


0.46 


48689 




Mil 


4.74 x 10 11 


110.3 


0.110 


0.106 


11.4 


542682 


4.16 x 10 10 


6.6 


1.29 


289205 


5.2 




6.99 x 10 10 


58.3 


0.071 


0.035 


9.94 


80000 


2.88 x 10 9 


1.6 


0.19 


20000 




M12 


5.65 x 10 10 


54.3 


0.076 


0.103 


7.21 


60000 


4.65 x 10 9 


2.6 


0.35 


20000 


17.7 




3.11 x 10 11 


95.9 


0.056 


0.075 


7.41 


330909 


1.96 x 10 10 


3.1 


0.46 


84596 





Another observational procedure to determine 
the structural parameters is by fitting the growth- 
curve of the luminous component (e.g., Burstein 
et al. 1987, Prugniel & Simicn 1997, Binggeli & Jcr- 
jen 1998). 

The accumulated luminous mass (or growth- 
curve) for a Sersic law is 



1999). The total luminous mass is given by 
27m 



Mr 



b 2, 



-r(2n) Eq-R, 



2 

e ' 



(7) 



M L (R) = 2n J Z(R)RdR 



^S i? c 2 7 (2n,6x 1 A 



(6) 



where x = R/R e , and (for a > 0) 

l(a,x) = [ 
Jo 



is the incomplete gamma function (Ciotti & Bcrtin 



were T is the complete gamma function. We as- 
sume here the approximation b(n) = 2n — 1/3 + 
4/(405n) + 46/(25515n 2 ), that provides a relative er- 
ror of < 10~ 7 in the range n G (1, 10). 

The structural parameters are somewhat depen- 
dent on whether a profile or growth-curve is fitted. 
Hence, we have fitted both a density profile, S(i?), 
and a growth-curve, M\,(R), to the luminous com- 
ponent of our remnants. We obtained the param- 
eters R e and Ml by minimizing the x 2 using the 
Levenberg-Marquardt method (Press et al. 1992). 

The central velocity dispersion of all the luminous 
particles was computed inside a circular region of 
projected radius R c /8; a standard practice in obser- 
vational studies (e.g., J0rgensen, Franx, Kjaergaard 
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0.1 0.5 1 5 10 



R [kpc] 

Fig. 1. Surface density profiles in "magnitudes" for all 
of our merger remnants; normalized to their "central" 
value. The outer limit of the fitting interval is taken 
to be 3-Ri, where R\ is the actual half- luminous mass 
radius. 

1996). Thus, the particular value of oo depends on 
the R c obtained from the adopted method of fitting. 
We computed ctq as 



ff0 = \ /V3l5>^~ <y * ))2; (8) 
\ »=i 

where N corresponds to the number of particles in- 
side the aperture, V Zi is the line-of-sight velocity of 
the z-th particle, and (V z ) is the mean velocity inte- 
grated along the line-of-sight. 

In order to have a better statistics and mimic 
observations, we looked at each remnant along 100 
random different lines-of-sight. For each projection 
we computed R c , Ml and cto as described above. 
This conforms our data set over which the linear fits 
are used to determine the scaling indexes. 

3.3. Interval of Fitting 

The structural parameters depend on the re- 
gion where the fit of the light profile is done (e.g., 
Caon et al. 1993; Graham 1998; Bertin, Ciotti & 
Del Principe 2002). Observationally, the significant 
radial range for the fit goes, for example, beyond the 
region dominated by seeing effects to that where the 
data is considered reliable. 

We have considered two inner radii £ for the fit- 
ting interval. First, we take a value of £ = 100 pc s» 
3e<j; the numerical resolution of our simulations for 



the luminous component. Secondly, we recalculated 
the scalings now adopting a value of £ = 1 kpc. As 
a reference, in a flat universe (Sl m = 0.3, Q\ = 0.7) 
with Hubble's parameter h = 0.7 an angular size of 
1".5 at, for example, the Coma Cluster [z — 0.023) 
corresponds to a physical size of ~ 700 pc. 

Contrary to observations, in iV-body simulations 
we can sample the complete luminous component 
and hence determine accurately the half-light ra- 
dius R\ of the remnant; that would correspond under 
ideal circumstances to the definition of R . However, 
given that R c is obtained from a fitting procedure, 
its value will be dependent on the radial range cov- 
ered by the fitting. In order to minimize such bias, 
we have done fittings inside a circular aperture of ra- 
dius 3R\. At this outer radius the range of sampled 
"magnitudes" [fi = — 2.51og£(i?)] of the luminous 
component is in average ss7; see Figure ^ This in- 
terval in magnitudes is similar to the found in some 
observational studies of ellipticals (e.g., Prugniel & 
Simien 1997, Bernardi et al. 2003a). 

3.4. Fitting Procedure 

To compute {a,f3,j} in the expressions of a 
least-squares linear fit, in log-space, was performed. 
Different methods exist to do such fit, and it is known 
that the slope of the fitted line depend on the fitting 
procedure (e.g., Feigelson & Babu 1992). 

Here the ordinary least-squares (ols) and the or- 
thogonal distance regression (odr) procedures are 
used. The ODR method is particularly useful when 
there is no clear distinction between the dependent 
or independent variable, and when both variables 
contain uncertainties. Moreover, the ODR fit is in- 
sensitive to whether the data are weighted or not 
(Wu, Fang, Xu 1998). The ODR line fitting is de- 
scribed in Appendix A. 

In order to estimate the standard deviations in 
the scaling indexes {a, (3, 7} a bootstrap technique 
(Efron & Tibshirani 1993) was used on our data set. 

4. SCALING RELATIONS 

In this section, we present the scaling relations, 
among the physical quantities that appear in the 
virial relation and their corresponding observational 
counterparts, obtained from our merger remnants. 

4.1. Velocities 

In Figure El we plot, in log-space, V 2 against a\ 
for each projection (dots) of the remnants, where 
the average value of CTg is indicated by a larger sym- 
bol. An homologous behavior for the velocity scaling 
would require that V 2 cx CTq. 
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Fig. 2. Virial velocity against the central velocity dis- 
persion for each projection (dots) of a remnant; aver- 
age values are indicated by solid squares. The solid line 
corresponds to an OLS fit to the data obtained using a 
Sersic growth-curve technique and £ = 100 pc; the scal- 
ing expression for this case is indicated. The broken line 
represents the fit by ODR. The slope of each line is indi- 
cated. 

A relation of the form 

V 2 cx a% , (9) 

in log-space, was done to out data set. In Table 2 
we list the results obtained for the a index under 
both fitting procedures ODR and OLS used, and both 
inner boundary radii £. The error in a by the OLS 
procedure is ±0.01, while that for the ODR procedure 
is ±0.02. 

It can be noted that both fitting procedures (OLS 
and ODR) lead to different values for the a index. 
The ODR procedure gives a value of a about 50% 
higher than that obtained by OLS, under all the con- 
ditions used for the fittings. 

Increasing £ in the fit leads to a small change 
(< 5%) in the a index for the same fitting proce- 
dure. At a fixed £, using the different fitting func- 
tions (profile or growth-curve), results in small vari- 
ations (< 5%) in the value of a. 

Averaging the results in Table 2, over £ and type 
of profile considered, we have for the OLS procedure 
a (o)ois = 1.51 while the ODR approaches a value 
of (a)odr = 2.35. In general, the values of a in our 
remnants deviate by ~ 20% from the homology ex- 
pected value of a = 2. It follows from these results 
that our merger remnants do not satisfy kinematical 



TABLE 2 
VELOCITY SCALINGS V 2 cx a'. 





-ftp 


pl/4 
n CG 


s P 


ScG 


FIT 


100 pc 


1.53 
2.40 


1.58 
2.40 


1.52 
2.36 


1.52 
2.35 


OLS 
ODR 


1000 pc 


1.46 
2.30 


1.54 
2.40 


1.48 
2.33 


1.46 
2.30 


OLS 
ODR 



i? 1 / 4 and S indicate a de Vaucouleurs and Sersic profile, re- 
spectively. Subscripts P and GC refer to whether a fitting to 
the surface density profile E(-R) or the curve of growth M^(R) 
is done, and OLS and ODR refer to the type of least-squares fit- 
ting procedure used. 



homology. 

We recall that the total (random plus rotational) 
kinetic energy of the system is related to V 2 , and 
ctq is essentially a measure of random motion. The 
rather small departure from homology of a seems to 
suggest that the contribution from rotational energy 
might be small in the luminous part of the remnants. 
This result appears consistent with some observa- 
tions of the rotational contribution to kinetic energy 
in ellipticals (e.g., Prugniel & Simien 1994). 

However, to test the degree of rotational en- 
ergy in our remnants a detailed kinematical analysis 
would be required, a topic that is under investigation 
at this time and to be presented in a future work. 

4.2. Radii 

In Figure |21 we plot the gravitational radius, i? g , 
of each merger remnant versus their corresponding 
effective radii, R Cl obtained from a fit using their 
Sersic growth-curve with £ = 100 pc, and considering 
all the randomly generated projections in the sky. 

A fit in log-space of the form 

R g cx Rg (10) 

was done. In Table 3 we list the different values of (3 
found under the different fittings conditions. Errors 
in j3, under both the OLS and ODR, are ±0.01. 

The (3 values do not change (~ 20%) as much 
as when the fit is done with an OLS or an ODR, in 
comparison to the « 50% change of the a index. An 
increase in (3 of w 20% is obtained when the fitting 
is done on the profile using a £ = 1 kpc instead of 
100 pc. We note that when the growth-curve method 
is used the {3 values decrease. 

The (3 index strongly depends on the adopted 
form of the fitting law, irrespectively of the fitting 
procedure. Its variation can be as much as w 200% 
when a Sersic law is used instead of the i? 1 / 4 law. 
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Fig. 3. Similar as to Figure 2, but for the gravitational 
radius R g and the effective radius Re. 

TABLE 3 

RADIAL SCALINGS R G cx 





rip 


pl/4 
n CG 


s P 


<ScG 


FIT 


100 pc 


0.09 


0.16 


0.27 


0.27 


OLS 




0.11 


0.19 


0.33 


0.33 


ODR 


1000 pc 


0.12 


0.11 


0.23 


0.23 


OLS 




0.13 


0.12 


0.26 


0.25 


ODR 



This behavior is most probably related to the one 
observed in observational studies (e.g., Khosroshashi 
et al. 2004), where totally different values of R e can 
be obtained if the brightness profile is not fitted by 
a suitable model. Here, the election of the adopted 
profile is reflected on the values of the /? index. 

Averaging all the values listed in Table 3 related 
to the de Vaucouleurs law, irrespective of the fitting 
procedure or £ value, we obtain a mean of (/3)dcV — 
0.13, while for the Sersic law a value of (f3)s — 0.27 
is obtained. 

The dispersion of values, due to projection ef- 
fects, around the mean value can be appreciated in 
Figure|3 where a Sersic growth-curve method is used 
to determine R c . Although not shown, the disper- 
sion around the mean value is somewhat larger when 
a i? 1 / 4 profile is fitted. 

The results obtained for the j3 index, under all 
the fitting conditions considered, indicate a strong 
breaking of an homologous scaling between R g and 
Re. In average, the (3 index deviates » 80% away 
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Fig. 4. Similar to Figure 2, but for the total mass M and 
the luminous mass Ml inferred from the fit. 

from the homology value of /3 = 1. 

4.3. Masses 

In Figure^lwe plot the total mass of the remnant 
M versus its luminous mass Ml, where the latter 
was obtained using a fit to a Sersic growth-curve and 
£ = 100 pc. 

A relation of the form 

MocMl (11) 

is assumed for the fitting in log-space of the data 
set. In Table 4 we list the values of the mass scaling 
index 7 for all the different fittings considered here. 
Errors in 7 in both procedures are ±0.01. 

The 7 index shows a smaller dispersion against 
projection effects in comparison to the a or fj in- 
dexes. This is probably related to the fact that Ml 
is an integrated quantity. 

The mass index 7 shows the same tendency to 
increase its value when an ODR fitting procedure is 
used in comparison to the OLS one. Increasing the £ 
value has the effect of lowering the value of 7 in all 
cases considered, and is rather stable under the two 
laws of luminous matter distribution assumed (de 
Vaucouleurs or Sersic). Averaging all values under 
the ODR and OLS procedures we find a (7) dr = 0.76 
and a (7) is = 0.69, respectively. 

It is difficult to adequately transform our Ml 
to a luminosity L in order to compare our results 
with observations, this due to the uncertainties on 
the mass-to-light ratios of ellipticals and the dissi- 
pationlcss nature of our simulations. Nonetheless, if 
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TABLE 4 
MASS SCALINGS MaM, 





R 1/4 


i? 1/4 


s P 


ScG 


FIT 


100 pc 


0.69 


0.76 


0.72 


0.73 


OLS 




0.78 


0.84 


0.79 


0.79 


ODR 


1000 pc 


0.64 


0.62 


0.68 


0.67 


OLS 




0.72 


0.70 


0.74 


0.72 


ODR 



we consider that M oc M^' 8 , as provided by an ODR 
fitting procedure using a Sersic law with £ = 100 pc, 
then the following total mass-to-luminous mass ratio 
is obtained: 

M M L 8 , .-0 2 

— oc — ex M, °- 2 . 
M L M l l 

Assuming a constant mass-to-light ratio, the pre- 
vious result would imply that M/L oc L~ 2 . Ob- 
servationally it has been found, using a Sersic pro- 
file (Trujillo, Burkert & Bell 2004), that M/L oc 
£,°- 06 ± - 04 . When comparing this result with the one 
obtained for our merger remnants, it follows that a 
constant mass-to-light ratio is not appropriate to re- 
produce the observational results. 

Nevertheless, a constant mass-to-light ratio does 
not seem so unrealistic especially when mass-to-light 
ratios oc L~ 0A have been found for dwarf ellipti- 
cals (Peterson & Caldwell 1993). Such dependency 
is more closely reproduced, under a constant M/L 
ratio, using other values of Table 4. However, a con- 
sistent comparison with observations requires addi- 
tional physics not considered in this study. 

5. DISCUSSION AND CONCLUSIONS 

We have carried out twelve iV-body simulations 
of binary mergers of disk galaxies, constructed us- 
ing the galaxy formation model of MMW and with 
properties consistent with a Tully-Fisher realization 
at z = 1. For the merger remnants, scaling rela- 
tions among the physical quantities that appear in 
the virial theorem and their observationally related 
ones were obtained. 

In particular we looked for relations of the form: 
V oc crj, R g oc RP, and M oc M£. It is found 
that the scaling indexes are sensitive to the fitting 
procedure (odr or OLS), to the inner starting radius 
of the fitting region, £, to the kind of law assumed 
to follow the luminous matter (de Vaucouleurs or 
Sersic), and to whether a profile or growth-curve is 
used. The 7 index results to be the more stable under 
all these different fitting conditions. 



In general, our results show that a strong break- 
ing of homology occurs in dissipationlcss mergers. 
We find that the a and 7 indexes are more sensitive 
to the fitting procedure, obtaining for the OLS pro- 
cedure a (a)ois = 1-51 and a (7) is = 0.69, while for 
the ODR procedure (a) dr = 2.35 and (7)odr = 0.76. 
The (3 index results to be more sensitive on the as- 
sumed law for luminous matter distribution, values 
(/3)dcv = 0.13 and (f3) s = 0.27 are found. 

An immediate consequence of our results is the 
existence of a non-linear scaling between the virial 
theorem, M oc V 2 R g , and its observational anal- 
ogy, Ml oc <JqR c . Moreover, the "constant" of pro- 
portionality between the physical and observational 
virial relations depends on the fitting procedure and 
the radial range where the fit is done. This indi- 
cates that care has to be taken when trying to obtain 
physical information from these observational pa- 
rameters; for example, in determining the mass-to- 
light ratio of ellipticals using a kinematical approach 
based on the observational virial relation (e. g., Pad- 
manabhan et al. 2004). Our results suggest that a 
dynamical mass estimate based on an homologous 
relation of the form Md yn oc a 2 R c is likely to be 
incorrect. 

Properties of remnants depend on the angular 
momentum and energy of the orbit of the progenitors 
(e.g., Naab & Burkert 2003, Gonzalez-Garci'a & Bal- 
cells 2005, Boylan-Kolchin, Ma & Quataert 2005). 
We expect that the scaling indexes, that reflect the 
breaking of homology in remnants, will also depend 
on these quantities. Boylan-Kolchin et al. (2005) ad- 
dress in an approximate manner, and with iV-body 
simulations, the degree of homology breaking and 
find a dependence on the type of orbit considered. 
Although they just used one pericentric distance and 
a radial orbit, aside of using only spherical mod- 
els, their result hints toward a dependency on an- 
gular momentum. Considering the small parameter 
space of simulations sampled out here, the values of 
{a, (3, 7} obtained in this study are to be taken as in- 
dicative of the values they might attain under more 
general conditions. A complete study will require 
to study the dependence of the scaling indexes on 
the energy and angular momentum of the orbit, and 
needs to be addressed in the future. 

Comparing our iV-body results directly to obser- 
vations is restricted by the uncertainties in trans- 
forming the luminous mass, Ml, to luminosities, L. 
This is not an easy problem and would require to 
include, among other things, gas and stellar popula- 
tions evolution models in the simulations. It is not 
clear at this stage how the scaling indexes would be 
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Fig. 5. Schematic diagram indicating the quantities used 
for fitting a line to a set of data points by orthogonal 
distance regression (odr). 

affected by the inclusion of this new physics into the 
problem. However, in conclusion, it is clear from a 
purely iV-body point of view that homology is not 
satisfied in merger remnants of spiral galaxies. 

Acknowledgments 

This research was funded by CONACyT-Mexico 
Project 37506-E. An anonymous referee is thanked 
for important comments that helped to improve the 
presentation and content of this work. 

APPENDIX A 

Orthogonal Fit to a Line 

In the ordinary least-squares (ols) method errors 
in the "independent" variable are minimized. How- 
ever, if there is no clear distinction which variable 
is dependent or independent a natural choice is to 
minimize errors in the normal direction to the sur- 
face fitted. This is the idea behind the orthogonal 
distance regression fitting (odr). We describe here 
the procedure used to fit a line by this method. 

Let the line equation be described by 



L = Xn 



(12) 



where D is a unit length vector along the line and t 
some scalar. If Xj is the 2-dimensional vector of the 
n data points to be fitted, it can be written as 



X< = X + C*D + £iD J 



(13) 



where T) ± is a unit vector perpendicular to D and 
Ci, Si are some scalars; see Figure El 

Let Yj = Xj — Xo, then we have that 

Yi - CiD = . 



The function to be minimized in ODR is 

n n 

£(X , dh 5>? = ]T(Y, - oB) 2 

i=l i=l 

We can write this as 

n 

S(X ,D) = ^[Yf(l-DD T )Y l ] 

i=l 

The best fit Xo is obtained by setting 
dE 



(14) 



(15) 



dX 



= -2(l-DD T )^Y, = 0, 



that is obtained by setting J^ILi (Xj— Xo) = 0, hence 



1 - 
Xo = — / Xj 



(y) 



(16) 



To obtain the best fit D we find de minimum of 
the following equivalent expression to (|15fl 

£(X , D) = D T jE[( Y * ' Y ^ ~ Y,Yf]| D , 

and 1 is the unit matrix. The term in brackets is a 
matrix, so we write this expression as 

£(X ,D) = D T M(X )D. 



Writing explicitly matrix M(Xq) we have 
M(X ) = 



-5 2 



-8 2 

xy 

5L 



where the S xy element is given by 

n 

^v = E(^-^»(»i-<w» 



(17) 



(18) 



(19) 



From linear algebra we recognize the expression 
in l|17() as a quadratic form. Its minimum value 
is provided by the eigenvector corresponding to the 
lowest eigenvalue of matrix (|18|l . In particular, nor- 
malizing this eigenvector to unity leads us directly 
the components D = (p, q) , and hence the best fit- 
ted line l|12|) . The line equation in parametric form 
can be written as 

x = x + tp , y = y +tq, 

and, eliminating the parameter t, we have the best 
line fitted by ODR as: 



y=-x + (y x ) ■ 

P P 



(20) 



Eigenvalues may be found by using Jacobi method 
(Press et al. 1992). 
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